Resistive and rectifying effects of pulling gold atoms at thiol-gold nano-contacts 
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We investigate, by means of first-principles calculations, structural and transport properties of 
junctions made of symmetric dithiolated molecules placed between Au electrodes. As the electrodes 
are pulled apart, we find that it becomes energetically favorable that Au atoms migrate to positions 
between the electrode surface and thiol terminations, with junction structures alternating between 
symmetric and asymmetric. As a result, the calculated IV curves alternate between rectifying 
and non-rectifying behaviors as the electrodes are pulled apart, which is consistent with recent 
experimental results. 

PACS numbers: 73.63.-b, 73.63.RT, 61.46.-w 
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Stimulated by recent experiments Q, |^ |1] , the use of 

organic molecules as basic components of nanostructured 
circuits has been envisaged as a way to the further mini- 
tuarization of electronic devices. In order to achieve real 
devices it is required a good understanding of the features 
affecting the conductance of a single molecule attached 
to electrodes. The atomic structure of the contact re- 
gion may strongly affect the conductance and the current- 
voltage {IV) characteristics of a molecular junction. For 
instance, the choice of different functional groups con- 
necting the molecular core to the electrodes has been 
shown to lead to very asymmetric IV curves '3| . Surpris- 
ingly, rectification can also occur in symmetric molecules: 
in some experiments the careful manipulation of the elec- 
trode distances in a symmetric molecule-gold junction 
can result in asymmetric characteristics 5, 6]. Since 
the organic molecule is symmetric, the probable reason 
for the asymmetric /F curves is the presence of asymmet- 
ric contacts. Experimentally, it has been found that the 
force required to break a molecular junction of a single 
octanedithiol molecule attached to gold electrodes is sim- 
ilar to the value of the force required to break a Au-Au 
bond in an atomic gold chain (1.5 - 1.8 nN) This fact 
strongly suggests that the dithiol molecule can pull Au 
atoms off the electrodes when a external force is applied. 
Molecular dynamics simulations also indicate that a thi- 
olate molecule can pull gold atoms off a stepped surface 
[Tof . Therefore, asymmetric contacts could be formed by 
stretching the molecular junctions which in turn could be 
the source of such asymmetries in the IV curves. 

In this paper we have investigated the effect of the 
strain on the atomic structure and on the electron trans- 
port properties of molecular junctions composed by sym- 
metric dithiolated molecules attached to Au(lll) elec- 
trodes. We find that the conductance is not very sensitive 
to small applied forces, in agreement with experimental 
observations made by Reichert et al. However, as the 
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FIG. 1: Alkanedithiol (Octanedithiol), Alkenedithiol 
(Octanedithiol), and [9,10-Bis{{2'-para-mercaptophenyl - 
ethinyl)-anthracene] (BPMEA) 



electrodes are pulled apart, we find that it becomes ener- 
getically favorable that Au atoms migrate to positions be- 
tween the electrode surface and thiol terminations, with 
junction structures alternating between symmetric and 
asymmetric. As a result, the calculated IV curves alter- 
nate between rectifying and non-rectifying behaviors as 
the electrodes are pulled apart, which is consistent with 
the results of Reichert et al. p ■ 

Our methodology is based on density functional theory 
(DFT) yjj] within the generalized gradient approxima- 
tion (GGA) |2| and norm-conserving Troullier-Martins 
pseudopotentials Ts] as implemented in the Siesta 
method jl^ . The solutions of the Kohn-Sham equa- 
tions are expanded as a linear combination of atomic 
pseudo-wavefunctions of finite range. We used a ba- 
sis set previously used to study methylthiol function- 
alized gold nanoclustersp^ ll^ IT^ . The conductance 
and the current are calculated within a nonequilibrium 
Green's-function formalism at the same DFT level using 
the TranSiesta method [l9j |. 

We have considered the three organic molecules shown 
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FIG. 2: Studied gold-octanedithiol-gold junctions. Similar 
junctions of alkenedithiol were also studied. 

in Fig. 1: an insulating alkane chain (octanedithiol); 
a conducting alkene chain (octenedithiol); and the 
symmetric conjugated polymer [9,10-Bis((2'-para- 
mercaptop/ienj/Z-ethinyl)-anthracene] (BPME A) . The 
molecules are connected to the gold electrodes through 
the sulphur atoms of the thiol groups, as shown in 
Fig. 2. The electrodes are described as in Ref. T^, with 
a (111) gold slab that connects to the bulk of gold leads 
through their self-energy matrix. The thickness of the 
slab is 9 layers, out of which one is defined as buffer 
[l9l |. and the lateral supercell repetitions used are 3x3 
and 5x5. We find that the 3x3 cells are large enough for 
the convergence of the total energy calculations, while 
the conductance calculations demand 5x5 cells for a 
convergence of the positions of the transmission peaks 
within 0.25 eVj2flj . 

In order to address the question of the stability of addi- 
tional gold atoms between the sulphur atom and the gold 
surface, we considered the five systems shown in panels 
(i) through (v) in Fig. [3 In (i) the molecule (alkane in 
the figure) is directly connected to the gold slabs; in (ii) 
an extra gold atom is placed in between one sulphur atom 
and the gold surface; in (Hi) extra gold atoms are found 
at both contacts; in (iv) two atoms are placed between 
the sulphur and the surface in one of the contacts; in (v) 
one contact has two extra atoms and the other contact 
has one. Similar configurations have been also considered 
for the alkene chain. The Helmholtz free energy is cal- 
culated (at T = 0) as a function of length and compared 
in Fig. 3 for the systems described, is the relevant 
energy for the experimental realizations of these devices 
since, even if the force is measured, it is the length what 
is defined by the experiment 9]. In our simulations, we 
increase the electrode-electrode distance in quasi-static 
steps by increasing the size of the unit cell in the z-axis. 
In each step the lattice vectors are kept fixed and the ge- 
ometry of the system is allowed to relax. The Helmholtz 
free energy is calculated according to 

^ — Etotal — ^J■AuNAu, (1) 



0( 



' 1 ' 1 ' 1 / 1 _ 
(a) / //N 

" / •■^iv' / 2.4 nN- 

- cTX^ 2.4 nN 1 
2.3 nN H 


' 1 ' 1 ' 1 ' 1 _ 

^ (b) / ^ 

/ 2.5 nN / "^-^Q^^ ' 

f 2.0 nN - 

Z j^^^^ 1.9 nN r 
1.9 nN 


1 1 1 1 1 1 1 A 

- (C) -v-^ ^ 

1 .^.8 nN/ 
/ el"*'/^ 2.6 nN- 

- i^^^ 2.2 nN 


1 1 1 1 1 1 1 1 , 

^ (d) / ^ 

/ 2.4nN/ * / 

1 / y .^/^ : 



L-Lq(A) 

FIG. 3: (Color online). Helmholtz free energies Q. of 
gold/molecule/gold systems as a function of the electrode- 
electrode distance. The top panels [(a),(b)] correspond to 
the gold/octanedithiol/gold systems, and the bottom panels 
[(c), (d)] to the gold/alkene/gold systems. In the left [(a),(c)] 
panels Q, was calculated with ^au ~ fJ-buik and the right pan- 
els [(b), (d)] with hau = fJiadatom- In each panel, the curves 
from left to right (circles — > squares — > triangles cross 
diamonds) correspond to the contact geometries from the left 
to the right in Fig. 2 [i ^ ii Hi ^ iv ^ v). 



where Etotal is the calculated total energy, Nau is the 
number of extra gold atoms, and iiau is the chemical 
potential of gold. As a lower bound for ^au we consider 
fJ-Au = l^-^Au^ total energy per atom for bulk gold. 
As a plausible upper bound we choose ^au = MaIT'"™' 
the energy of an adatom on the Au(lll) surface, i.e., 
the energy difference between the Au(lll) slab with and 
without an adatom. 

Fig. 3 displays the Helmholtz free energy as a function 
of the length. We also indicates the value of the force 
at the intersections (the force is Fz = —SU/Sz as coming 
from the low L side) . These force values can be seen as es- 
timates of the forces required to extract atoms from the 
gold surface, at least in the adiabatic low-temperature 
limit. Indeed, the force required to extract the first atom 
(1.8-2.2 nN) is iii good agreement to the value obtained 
by Kriiger et aL\W^ using Car-Parrinello molecular dy- 
namics simulations. These extracting forces are not very 
different for alkane and alkene junctions. This suggests 
that the forces required to pull off Au atoms from the sur- 
face do not depend much on the organic molecular core, 
and one could expect similar values for other dithiolated 
molecules. The structures (iii) and (iv) in Fig. |21 have 
the same number of atoms. Their Helmholtz free ener- 
gies shown in Fig. |21 indicate that structure (iii) is the 
most stable one, by 0.7 eV. Therefore, considering ener- 
getics aspects only, the results of Fig. 3 would indicate 
the following sequence of junctions during the pulling of 
the electrodes: (i), symmetric; (ii), asymmetric; (Hi), 
symmetric; (v), asymmetric. Since an additional gold 
atom between the sulphur and the electrodes generates 
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FIG. 4: (Color online). Conductance through octenedithiol 
as a function of the inter-electrode separation: black circles 
correspond to the Fig. 2{i) junction, in which the molecule 
is directly connected to the Au(lll) surfaces; blue squares 
are for the Fig. 2{ii) junction, in which an extra gold atom is 
placed between a sulphur atom and the electrode; red trian- 
gles are for the Fig. 2{i) junction, for unrelated electrodes. 



asymmetries in the IV characteristics of dithiol junc- 
tions, as we shall show below, the sequence of symmetric 
and asymmetric contacts described above is consistent 
with the observation of both symmetric and asymmetric 
IV curves by Reichert et al. ^ for a BPMEA attached to 
gold slabs. In that experiment, starting with a symmetric 
IV they increased the inter-electrode distance obtaining 
an asymmetric IV. After further manipulations they re- 
covered the symmetric IV followed again by asymmetric 
curves. 

Fig. 0] shows the conductance of the octenedithiol as a 
function of the distance between electrodes for the first 
and second junctions in Fig[21 i.e., with and without an 
extra gold atom (circles and squares, respectively). There 
is no appreciable change in conductance for any of the 
systems over the quite substantial extensions displayed, 
which is consistent with what is observed experimentally 
by Reichert et al. ^ . They find essentially unaltered IV 
characteristics upon extension of up to a few A, and sug- 
gest that the strain is relieved by the soft gold bonds. To 
test this hypothesis we performed calculations keeping 
the geometry of the gold electrodes fixed, relaxing only 
the molecule for each elongation. In this non-relaxed 
case, the conductance of the alkene chain increases from 
approximately 0.4 Gq to almost 1 Go when the electrode 
distance is increased by only 0.6 A, as shown Fig.^ This 
shows that the weak dependence of the conductance with 
elongation is due to the relaxation of the electrode geom- 
etry near the contact region. 

In contrast with the small effect of elongation on the 
conductance described above, the effect of extra Au 
atoms in the conduction of the molecular junction is very 
large, as indicated in Fig. 0| The figure shows that the 
conductance of the Au/octenedithiol/Au junction is re- 
duced by a factor of three upon the the inclusion of an 
extra Au atom between one of the molecular thiol termi- 
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FIG. 5: (Color online). Total transmission coefficient ver- 
sus energy, (a) and (b) are for BPMEA, (c) and (d) are for 
octenedithiol. The upper panels compare the junctions with- 
out (dashed lines) and with (continuous lines) an extra gold 
atom at a thiol-electrode junction. The lower panels show 
the effect of the applied bias on the extra-atom junctions in 
both cases. In lower panels, the energy origin corresponds to 
the average of the chemical potentials in the left and right 
electrodes, iii±ii«. 



nation and the electrode. This effect is further explored 
in Fig. 5 for BPMEA [(a) and (b)] and for octenedithiol 
[(c) and (d)]. The upper panels show the effect of the 
additional gold atom in the low bias conductance, while 
the lower panels show the asymmetric bias dependence 
for the junctions with the extra atom. The effect of the 
extra atom on the conductance is significant, decreasing 
G to approximately one third of its original value for the 
alkene, while it increases from 0.02 Go to 0.07 Go in the 
case of BPMEA. It is worth to note that the calculated 
values of conductance of the BPMEA junction without 
an additional gold atom(0.02Go) is in qualitative agree- 
ment with the experimental result obtained by Reichert 
et alQ, ~ O.OIGq. The apphed bias produces a shift 
of the spectral features in the transmission, which gives 
rise to asymmetry in the IV characteristics. Under ap- 
plied bias values of 0.1 V and -0.1 V, we obtain currents 
of 1.2 ^A and 0.7 ^A, respectively, for BPMEA. Fig. 
shows the rectification ratio, for the alkene junction 
with and without the extra atom, displaying a clear rec- 
tification in the former case (circles). The much smaller 
rectification observed in Fig. 6 for the case without extra 
Au atoms is a result of the absence of a left-right symme- 
try in the device (the molecule has inversion symmetry 
while the contacts have reflection symmetry), which gives 
rise to a small intrinsic effect. The introduction of the 
extra atom produces an increase of an order of magnitude 
in the rectification. 

In conclusion, we find that the forces normally present 
in experiments with molecular nanojunctions do not 
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FIG. 6: (Color online). Rectification ratio ^ for the 
octenedithiol junctions in Fig. 2(i) (direct contact to the elec- 
trodes, black squares) and Fig. 2{ii) (one extra Au atom in 
one of the thiol junctions, red circles). 



directly affect the conductance characteristics through 
these devices due to the softness of the gold leads. They 
can, however, stabilize extra gold atoms at the electrode- 
thiol junctions, which do change the conductance and 
give rise to asymmetry and appreciable rectification ef- 
fects, which is in qualitative agreement with recent ex- 
perimental results. Our results concerning the mechani- 
cal properties of the junctions are weakly dependent on 
the actual molecule between the electrodes, being mainly 
determined by the thiol-gold chemistry. 
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